Three-dimensional topological phase on the diamond lattice 
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An interacting bosonic model ol Kitaev type is proposed on the three-dimensional diamond lattice. 
Similarly to the two-dimensional Kitaev model on the honeycomb lattice which exhibits both Abelian 
and non- Abelian phases, the model has two ("weak" and "strong" pairing) phases. In the weak 
pairing phase, the auxiliary Majorana hopping problem is in a topological superconducting phase 
characterized by a non-zero winding number introduced in A. P. Schnyder, S. Ryu, A. Furusaki, 
and A. W. W. Ludwig, Phys. Rev. B 78, 195125 (2008) for the ensemble of Hamihonians with both 
particle-hole and time-reversal symmetries. The topological character of the weak pairing phase is 
protected by a discrete symmetry. 



I. INTRODUCTION 

The recent discovery of Z2 topological insulators, a 
band insulator with particular topologicaliJjjhajcact^ of 
Bloch wavefunctions, came as a surprise. BHtJo IjIjiJa'El On 
the one hand, Z2 topological insulators are close relajtiypa 
to more familiar integer quantum Hall (IQH) states £3'Lll 
As in an IQH state in the bulk, they are characterized by 
a topological invariant (Z2 invariant). As in an IQH state 
with boundaries, they support stable gapless boundary 
states that are robust against perturbations. On the 
other hand, unlike the IQH states, time-reversal symme- 
try (TRS) is a prerequisite to the existence of Z2 topo- 
logical insulators. In fact, as soon as the TRS of a Z2 
topological insulator is broken, it becomes possible to 
deform in a continuous manner a band insulator with a 
trivial Z2 topological number into one with a non-trivial 
Z2 number. 

Time- reversal symmetry for spin 1/2 particles is not 
the only discrete symmetry for which a topological dis- 
tinction of quantum ground states arises. A system- 
atic and exhaustive classification of topological band in- 
sulators and mean-field superconductors has been pro- 
posed in Ref. |lj by relying on the discrete syirmetries 
of relevance to the theory of random matrices.EHlEj In 
three spatial dimensions, it was shown that, besides 
the Z2 topological insulator in the symplectic symmetry 
class, there are precisely four more symmetry classes in 
which topological insulators and/or superconductors are 
possible.E§E3 For three out of the five symmetry classes 
of random matrix theory, we introduced a topological in- 
variant ly (winding number), which distinguishes several 
different topological insulators/superconductors, just like 
the Chern inteffiarj-distinguishes different IQH states in 
two dimensions il3'EJ 

While the classification given in Ref. [l| is for non- 
interacting fermionic systems, strong correlations among 
electrons (or spins) might spontaneously give rise to 
these topological phases, by forming a nontrivial band 
structure for some, ppssibly emergent, fermionic excita- 
tions (e.g., spinous). til It is the purpose of this paper 
to demonstrate how topological insulators (superconduc- 
tors) emerge as a result of strong correlations. We will 



show that it is possible to design an interacting bosonic 
model with emergent Majorana fermion excitations, the 
ground state of which is a topological insulator (super- 
conductor) with non-vanishing winding number, ^ 0. 

Our model is a natural generalization of the spin- 1/2 
model on the honeycomb lattice introduced by Kitaevta 
to the three-dimensional diamond lattice with four di- 
mensional Hilbert space per site. The Kitaev model on 
the honeycomb lattice has two types of phases: the so- 
called Abelian and non- Abelian phases, r^he Abelian 
phase is equivalent to the toric code modeO and an ex- 
actly solvable model proposed by Wen,E2l which in turn 
is described by a Z2 gauge theory. On the other hand, 
the non- Abelian phase is in the universality class of the 
Moore-Read Pfaffian state. Each phase corresponds to 
the weak and strong pairing phases of two-dimen^nal 
spinless chiral p-wave superconductor, respectivcly,Eil the 
latter of which is an example of a topological supercon- 
ductor in symmetry class p .of. Altland-Zirnbauer classi- 
fication in two dimensions Otj'Efl 

Similarly to the Kitaev model on the honeycomb lat- 
tice, the ground state of our model can be obtained from 
a Majorana fermionic ground state (with a suitable pro- 
jection procedure). Our model has two phases, which 
we also call strong and weak pairing phases. In par- 
ticular, in the weak pairing phase, the ground state is 
given by a topological superconducting state in symme- 
try class Dili of Altland-Zirnbauer classification, and in 
the universality class of a three-dimensional analogue of 
the Moore- Read Pfaffian state discussed in Ref. T 
B phase of "^He is also in this universality classJ 
The topological character of the ground state is protected 
by a discrete symmetry transformation, which is a combi- 
nation of time-reversal and a four-fold discrete rotation, 
the latter of which forms a subgroup of a continuous U(l) 
symmetry of our model. Spin-1/2 models of Kitaev type 
on the diamond lattice, and OH^- nthpr , three-dimensional 
lattices have been constructed HJ'OES For these models, 
however, there is no phase analogous to the non- Abelian 
phase in the original Kitaev model, and the ground states 
discussed there have a vanishing winding number. Ex- 
tensions of the spin 1/2 Kitaev model to models with 
4-dimensional Hilbert space per site have been studied 
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in Rcfs. |20|j2^,||,|9|j3^. 



II. LOCAL HILBERT SPACE AND DISCRETE 
SYMMETRIES 

We start by describing the local Hilbert space of our 
model, defined as it is at each site of some lattice. Con- 
sider the four-dimensional Hilbert space spanned by the 
orthonormal basis 



This is nothing but the charge conjugation matrix C = 
47^7'^ for the gamma matrices in the Dirac representa- 
tion. As 9 is TRS for half-integer spin, 6^ = -1. Note 
also that 



(7) 



Second, if the local Hilbert space is interpreted as de- 
scribing two spin 1/2 degrees of freedom, we can consider 
a TR operation T defined by 



cr = ±l, T = ±l. 



(1) 



This space can be viewed, if we wish, as describing the 
four-dimensional Hilbert space of a spin 3/2 degree of 
freedom, or as a direct product of two spin 1/2 Hilbert 
spaces. In the latter case, one can view these two spin 
1/2 degrees of freedom as, say, originating from spin and 
orbital. We will denote two sets of Pauli matrices, cr'^ = 
crO,cr^,cr^cr^ and T^' ^ (/i = 0,1, 2, 3), each 

acting on a and r indices, with and being 2x2 
unit matrices. 

We shall represent the Hamiltonian in terms of two 
sets of Dirac matrices (Dirac representation). 



= /3 = 7°, (2) 



and (^M-o, 1,2,3 (^gj^jj-al representation), 



(3) 



where a = 1, 2, 3. The two sets {a^^} and {C^} are related 
to each other by = ia^^i^^^^ = ia^{cr'^ (8) r**), and 
satisfy the Dirac algebra, 

{«^a''}=={c^r}-25^^ m,^^-o,...,3. (4) 



A. discrete symmetries 

In the following, wc will consider three antiunitary dis- 
crete symmetry operations, T, T', and 8. They are char- 
acterized by 



-1, e-^ 



-1. T" 



-1. 



(5) 



T = {ia«) ® {iT«)K, 

Ta^'T-^ = -cr°, Tt^'T-^ = -r°, (8) 

with a = 1, 2, 3. Note that = +1. Under T, a and C 
arc transformed as 



A" 



(9) 



where covariant and contravariant vectors are defined as 
at' = (/?, a") and = (/3, -a"). 

As we will see later, while the cr-part of our Hamil- 
tonian is fully anisotropic in a space, the r-part of the 
Hamiltonian is invariant under a rotation around t"^ axis. 
In particular, it is invariant under a rotation i? by 7r/2 
around r** axis. 



i? = , i? = (r" + ir'')/\/2. (10) 



Under i?, a and C arc transformed as 



Ri-f-PR-^ = +i7^7°. 



(11) 



By combining T with R we can define yet another antiu- 
nitary operation, T' = i?T, 



T' ^RT^ [iry - r°)iCT'^A7\/2, 



(12) 



In the sequel, we will treat two distinct operations for 
time-reversal (TR). 

First, if the local Hilbert space is interpreted as de- 
scribing a spin 3/2 particle, the natural TR operation G 

IS given by 6 = r^e-"^y K, where rj stands for an ar- 
bitrary phase (will be set to one henceforth), Sy^"^ is a 
four by four matrix representing the y component of spin 
with S = 3/2, and K implements the complex conjuga- 
tion, KiK~^ = —i. If wc take \aT) to be the basis that 
diagonalizes S"^ (magnetic basis, |3/2,to)), then 
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(6) 



Below, with a slight abuse of language, we will call this 
operation T' time-reversal operation (TR) . When applied 
to a and C,, 



(13) 



i.e., TRS T' exchanges a and (, and covariant and con- 
travariant vectors. Notice that 



T'2 = iT^, T" 



(14) 
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FIG. 1: (Color online) The diamond lattice and the six Ma- 
jorana fermions -'S Sites on the sublattice A (B) are 
denoted by an open (filled) circle. 



III. HAMILTONIAN 

In the two-dimensional Kitaev model, different types 
of interactions (represented by three 2x2 Pauli matri- 
ces) are assigned to three distinct types of bonds that arc 
defined by their different orientation in the honeycomb 
lattice. It is this anisotropic nature of the interaction 
that makes the Kitaev model exactly solvable. This con- 
struction can be extended to the diamond lattice where 
there are four distinct types of bonds with different ori- 
entation. 

The diamond lattice is bipartite and consists of two in- 
terpenetrating fee lattices shifted by a(— 1, 1, — 1)/4 along 
the body diagonal, where a is the lattice constant. We 
label sites rA and on the two sublattices A and B of 
the diamond lattice as 



rA 



i=l 



rs =rA + So, rui £ Z, 
where the primitive vectors are given by 



ai = 



02 




03 = 



(15) 



(16) 



We have also introduced the 3-componcnt vectors 

-1 

a \ \ a I 

Sl 



S3 




(17) 



which connect the nearest neighbor sites. (Fig. |i|). 
The Hamiltonian we study in this paper is defined by 

^ = -E^'' E ("X + cj'cn- (18) 

fi—O /i— links 



Here, the sites j and k are end points of a link of type fi. 
There are four types of links = 0, 1, 2, 3 in the diamond 
lattice which can be distinguished by their orientations. 
Hamiltonian ([l8| ) can also be written in terms of and 



as 



H 



fi—d — links 



CT,. CTj, [Tj Tj, 



(19) 



This Hamiltonian is invariant under discrete symmetries, 
T, R, T' and 0, and enjoys a U(l) symmetry for rotation 
around axis. 



IV. MAJORANA FERMION 
REPRESENTATION 

Let us consider the (local) Hilbert space in which we 
have spS|[Majorana fermions {A''}p=o 5 per site, which 
satisfyE2E3 



{AP,An = 2<5^'', p,g = 0,...,5. 



(20) 



To construct 4-dimensional Hilbert space out of the 8- 
dimcnsional Hilbert space, we introduce the fermion 
number operator by 



p=0 



(21) 



The eigenvalue of £> = ±1 can then be used to select a 
4-dimensional subspace of the full Hilbert space, which 
will be called the physical subspace. 

We can construct 15 generators of so(6) from the Ma- 
jorana fermions as 

TP" = iXP\\ p^q. (22) 

Within the physical subspace, can be expressed as 

a^=^^^ M = 0,1,2,3. (23) 

Since F^^ anti-commutes with a'', we make the identifi- 
cation 



The second set of the gamma matrices {C^} is 



(24) 
(25) 



The Majorana fermions naturally inherit the symme- 
try operations on a and C. The symmetry conditions are 
automatically satisfied if we define T,R, and T' opera- 
tions on Majorana fermions by 

TA^r-i = A^, TX''T-^ = X\ s = 4,5, 

i?A^i?-i = _A^, r(^Ar-^ = isy ( Jjs 



T'y^T'-^ = - A., T' f C.^ ) T'-^ = is'y 



A4 
A5 



A4 
A5 



(26) 
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Here, the definitions of the covariant (A'') and contravari- 
ant (A^) vectors follow from those of (C^) and 
(Cfi), and we have introduced another set of Pauli ma- 
trices sM=o, 1,2,3 acting on A^'^ with s° being 2x2 unit 
matrix. The discrete rotation operator R can be written 
in terms of the Majorana fcrmions A as e™^ A''/4^ 



V. SOLUTION THROUGH A MAJORANA 
HOPPING PROBLEM 

In terms of the Majorana fcrmions the Hamiltonian 
can be written as 



H 



/i — links 



{X^Xt + Af A|) . (27) 



where we have introduced a hnk operator by 

\ fj'jk (28) 



Ujk := iXj Xk ' 



with ^Ljk = 0,1,2,3,4 depending on the orientation of 
the link ending at sites j and k. Note that TR (T or T') 
operation flips the sign of a link operator. 



TujkT ^ = -Ujk, RujkR ^ = 
T'ujk.T'~^ = -Ujk. 



(29) 



(This is also the case for TRS on the link operators in 
the spin-1/2 honeycomb lattice Kitaev model.) When 
necessary, this sign flip can be removed by a subsequent 
gauge transformation for Majorana fcrmions A^'^ on ci- 
ther one of sublatticcs, if the underlying lattice structure 
is bipartite (see below). 

What is essential to observe is that all Uj^ appearing 
in the Hamiltonian commute with each other and with 
the Hamiltonian. They can thus be replaced by their 
eigenvalues Ujk = ±1, and the interacting Hamiltonian 
reduces to, for a fixed configuration of the Z2 gauge field 
{wj-fc}, a simple hopping model of Majorana fermions. 
Observe that both A''' and A^ Majorana fermions feel the 
same Z2 gauge field. The ground state of the Hamiltonian 
can then be obtained by first picking up the Z2 gauge field 
configuration that gives the lowest ground state energy 
for the Majorana hopping problem, and then projecting 
the resulting fermionic ground state onto the-|physical 
Hilbert space. According to Lieb's theorem,!^ the Z2 
gauge field configuration that gives the lowest ground 
state energy has zero Z2 vortex for all hexagons, and 
hence we can take Ujk = 1 for all links. 

For notational convenience, for a Majorana fermion at 
the jth site located at (rg) on the sublattice A {B), we 
denote af.^ (Kb ^j)- With periodic boundary 

condition and with the Fourier transformation, af. , = 



Efee"^"--4/^A^ and K 



where |Ayi_B| is the total number of sites on the sublattice 
A, i?, respectively, the Majorana hopping Hamiltonian in 



the momentum space is 



s k 

where we have defined 



n{k) 



-i$*(fc) 



M=0 



(31) 



and noted = aj, when fc 7^ 0. The energy spectrum 

E{k) is given by E{k) = ±yi$(W> with two-fold de- 
generate for each k. 



A. symmetries and topology of the Majorana 
hopping Hamiltonian 

We have reduced the interacting bosonic model to the 
Majorana hopping problem. This auxiliary Majorana 
hopping Hamiltonian is, in the termiftelegy of Altland 



and Zirnbauer, in symmetry class D)13'0 i.e., the en- 
semble of quadratic Hamiltonians describing Majorana 
fermions. (See Appendix ^.) In more general situations 
(which we will consider below), the auxiliary Majorana 
hopping Hamiltonian is given by 



k 



-fc' "—A;' -fe 



/at 



(32) 



where X describes a Hamiltonian for Majorana fcrmions, 
and satisfies 



X\k)^X{k), A''^(-fc) -A'(fc). 



(33) 



This is the defining property of symmetry class D. Be- 
low, to describe the 4x4 structure of the single-particle 
Majorana hopping Hamiltonian X{k), we introduce yet 
another set of Pauli matrices c^=0:^,'2-3 acting on sublat- 
tice indices. 

If our bosonic model further satisfies TRS T', the 
Hamiltonian X for the auxiliary Majorana hopping prob- 
lem respects 



(34) 



where the factor can be thought of as a gauge transfor- 
mation, adding a phase factor e™ for Majorana fermions 
on B sublattice 6*, and can be removed by a unitary 
transformation 6'^ — > —6''. With this further condition 
arising from T', the relevant Altland-Zirnbauer symme- 
try class is class DHL (See Appendix]^) 

In Rcf. ^ it has been shown that the space of all pos- 
sible quantum ground states in class Dili in three spa- 
tial dimensions is partitioned into different topological 
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sectors, each labeled by an integer topological invariant 
V. To uncover this topological structure and introduce 
the winding number, we observe that all Hamiltonians 
in symmetry class Dili can be brought into a block off- 
diagonal form. For X , this is done by a unitary transfor- 
mation 

[/^t/jt/i, X X ^ U^U^XUlul, (35) 
where the first unitary transformation rotates — > 

Uisyul = -s^ 

C/i = {s° - ^s")/V2, (36) 

whereas the second unitary transformation exchanges 
2nd and 4th entries, 

U2 = (s" + s")c"/2 + (s" - s--)c72. (37) 

The combination of Ui and U2 diagonalizes s^c^ as 

U^U^syc'Ulul = -diag (1, 1, -1, -1) . (38) 

After the unitary transformation, we find 

This block off-diagonal structure is inherited to the spec- 
tral projector P(fc), = P, which projects onto the 
space of filled Bloch states at each fc, 

2^-1 ^(^t(,) t'O' (^°) 

The integer topological invariant is therv-defined, from 
the off-diagonal block of the projector, astil 

'^M - / T^^'"' [{q-'d,q) ■ (q-'d^q) ■ {q-'d,q)] , 
JBz ^'^^ 

(41) 

where ^^v^p = kx,ky,kz, and, the integral extends over 
the first Brillouin zone (Bz).EHl The non-zero value of the 
winding number signals a non-trivial topological struc- 
ture, an observable consequence of which is the appear- 
ance of gapless surface Majorana fermion modes. 

VI. STRONG PAIRING PHASE 

When one of the coupling J^^ is strong enough com- 
pared to the others, the spectrum for the Majorana 
fermions is gapped. In this phase, the winding num- 
ber v is zero. This phase can be called "strong pairing 
phase" , following the similar phase in the BCS pairing 
model. Because of the trivial winding number, there is 
no surface stable fermion mode in the auxiliary Majo- 
rana hopping Hamiltonian, when it is terminated by a 
surface. The properties of this phase can be studied by 
taking the limit Jq 3> Ja=i,2,3 > 0, say, and then devel- 
oping a degenerate perturbation theory. In this limit of 



isolated links, there are four degenerate ground states for 
each link ending at the two sites (r^, r^) = {rA, rA + sq), 

^ (I T);j T);, + 1 K.KMrsK, m 

with ^ = ±1 where | • • • represents the state for 
(Tr, and Tr, respectively. The effective Hamiltonian acting 
on these degenerate ground states is defined on the cubic 
lattice, since each Jo link at {rA,rA + so) is connected 
to six neighboring links at ± a^, where Qi = sq + Si 
with i = 1,2,3, upto the 4th order in the degenerate 
perturbation theory. If we use notations, cr^ p^, — s- 
/i'', the effective Hamiltonian upto the 4th order in the 
degenerate perturbation theory is (upto constant terms) 

(y,z,x),(z,x,y) 

where p stands for a plaquette surrounded by r, r -f a, 
r + aj and r 4- -f- aj , and 

A similar model on the cubic lattice was discussed in 
Refs. |g|,|7lj2|. 

VII. WEAK PAIRING PHASE 

When all are equal, = J, the energy spectrum 
of the Majorana hopping Hamiltonian Eq. ( pO| ) has lines 
of zeros (line nodes) in momentum space. This gapless 
nature is, however, not stable against perturbations that 
respect TRS T'. This can be illustrated by taking a four 
"spin" perturbation in the gapless phase, defined on three 
sites j,k,l, where sites j and k are two different nearest 
neighbors of site I. Let us take, as an example, 

_ p04p24p35pl4 
j I I k 

^ iUjiXj X Di X uikXt, (45) 

where we take the link emanating from sites j (k) and I 
to be parallel to sq (si). If perturbations of this type are 
small enough, relative to the excitation energy of a Z2 
vortex loop (line), i.e., an excitation which flips the sign 
of Z2 flux threading hexagons, we can contain ourselves 
in the vortex- free sector where Ujk = +1- Thus, the 
above four "spin" perturbation leads to a next nearest 
neighbor hopping terms of the Majorana fermions. To 
respect TRS T', Eq. (E5|) can be supplemented with its 
TRS partner, C°afCfCr= T'a"j(fafalT'-\ leading to a 
perturbation 

HLn = E ^^'o H^f^ic^k - (« 0] 

= I E K^flkUjiuik (A^At. - \]Xl) , (46) 
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where the summation extends over all sites labeled by I 
and their nearest neighbors j and k, with the link em- 
anating from sites j (fc) and I parallel to (s^), and 
Kjif, £ R. Similarly, the following perturbation defined 
on three sites j,l,k is also allowed by TRS T', 

H„\„ = K^ik K(a^*7^°r)/C.^ + (« - 0] 

= I J2 K^ikumk (A|A| + X]Xt) (47) 

with K^n- G K- We choose K^i^ in such a way that these 
perturbations lead to the following next nearest neighbor 
terms in the Majorana hopping Hamiltonian (see Fig. ||) 



nnn ^ 

2i 



A„ s A 



TTX ^_\^\T T 
nnn rj ■ / j 



r+si -S3 



h.c, 

+ Ar + S2— S3 



(48) 



where = (A^, A^) and i^^-^ G R. 

With these perturbations, the Majorana hopping 
Hamiltonian in momentum space is given by 



X{k) = 



e(fc) i$(fc) 
-i$*(fc) -e(fc) 



(49) 



where <I>(fc) comes from the nearest neighbor hopping 
(p7|), whereas the off-diagonal part 0(fc) comes from the 
next nearest neighbor hopping terms ( ffS]) and is given by 

e(fc) = e"(fc)s" + e"(fc)s^ 



e^(fc) = K"" 



e^(fc) = iC^sin- 



(50) 



Observe that this Hamiltonian indeed satisfies class DHI 
conditions ( |3^ ) and (^4|). ft can then be made block 
off-diagonal by the unitary transformation J7, with the 
off-diagonal block being given by 

D{k) = -Im$(fc)s" + e^(fc)is^ 

+ Q''{k)isy + Ke^{k)is''. (51) 

The energy spectrum at momentum k is given by E{k) ~ 

With these perturbations, we can indeed lift the de- 
generacy except at three points k = Qx.y,z in the Bz, 
where 






-0.4 -0.2 0.2 0.4 



-0.4 -0.2 0.2 0.4 



FIG. 2: (Color online) (Top left) The choice of the second 
nearest neighbor couplings K^ik (blue links) and K^i^ (red 
links). (Top right) The phase diagram in term of A"^ and 5J\ 
with J2,3,4 = J = 2, and = 1. (Bottom) The numerical 
evaluation of the winding number as a function of the second 
neighbor coupling K'' and the distortion 5 Ji . In the left panel, 
the winding number is computed for K'^ = ±1 with changing 
SJi continuously, whereas in the right panel 5Ji is fixed, SJi = 
1. For if^ = 1 (A'^ = -1), v^l{v = -l) when (5Ji > 0. 



The dispersion around these points is Dirac-like, 

X{Qa + q)- Jqac" + K-{qb - qc)s''c' 
-{qy + qz)s^c\ 



+ - 



(53) 



(52) 



where (a,6, c) is a cyclic permutation of {x,y,z). These 
three-dimensional Dirac fermions can be made massive 
by further adding a slight distortion in the nearest neigh- 
bor hopping, Ji ^ Ji + (5Ji, say. This gives rise to, at 
Qx,y.z^ a perturbation to X{k) which takes the form of a 
mass term to Dirac fermions, —SJid'. 

For definiteness, we now set J^=2.3,4 = = 2, = 1, 
Ji = J + (5Ji, and vary and 5Ji. In the (JJi,if^)- 
plane, there are phase boundaries represented by 5Ji = 
and the half-hne K'' = with 5Jx > (Fig. |). On the 
line SJi = 0, the spectrum is Dirac like except at the ori- 
gin {6Ji,K^) = (0, 0) where the band gap closes at Qx,y,z 
quadratically in one direction in momentum space (a sim- 
ilar gapless point is discussed in Ref. |3^). To determine 
the topological nature of the three gapped phases, the 
first and second quadrants in ((5Ji, X^)-plane, and the 
region 5Ji < 0, we computed the winding number, by nu- 
merically integrating the formula Eq. (^). Integral ( |4l| ) 
quickly converges, to a quantized value = 0, ±1 as we 
increase the number of mesh in momentum space. While 
the winding number is identically zero when SJi < 0, it 
takes either = +1 or = — 1 in the phases SJi > 0, 
depending on the sign of K^. The complete structure 
of the phase diagram including the value of the wind- 
ing number is presented in Fig. ^j. In the phases with 
non-zero winding number, there appears a gapless and 
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stable surface Majorana fermion mode, when the Majo- 
rana hopping Hamiltonian is truncated by a boundary, 
signahng non-trivial topological character in the bulk. 



VIII. DISCUSSION 

We have constructed a three-dimensional interacting 
bosonic model which exhibits a topological band struc- 
ture for emergent Majorana fermions. We thus take a 
first step to explore topological superconductors arising 
from interactions, rather than given by some external 
parameters at the single particle level, such as external 
magnetic field or spin orbit coupling. Although the Ki- 
taev model does not look particularly realistic as it is 
anisotropic both in real and spin spaces, it has played an 
important role in deepening our understanding on two- 
dimensional topological order. (See, for example, Refs. 
||,||j35|j3^,|^,|||3|,|^). Also, there have been a pro- 
posal to realize the Kitaev m odel , in terms of cold po- 
lar molecules on optical Iatticegl3e3 and superconducting 
quantum circuits O Interactions which are anisotropic 
both in real and internal spaces can appear in systems 
with orbital degrees of freedom, such as the orbital com- 
pass model. Indeed, it is worth emphasizing that our 
model, in the absence of four spin interactions, pos- 
sesses a U(l) rotation symmetry unlike the original Ki- 
taev model and its variants. Thus, say, identifying r as 
a spin 1/2 degree of freedom and a as an orbital degree 
of freedom, it might be realized as a XY analogue of the 
Kugel-Khomskii model. Finally, while our model is de- 
signed to have a Gutzwiller-type projected wavefunction 
as its exact ground state, such ground state wavefunc- 
tions can appear in much wider context, which can be 
explored, e.g., in terms of a variational approach with 
slave particle mean field theories .CJ 
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APPENDIX A: CLASS Dili SYMMETRY CLASS 

In this appendix, we review the symmetry classification 
of the Bogoliubov-deXfennes (BdG) Hamiltonians by Al- 
tland and ZirnbauerJ13 which is relevant to our auxiliary 
Majorana hopping problems. We consider the following 
general form of a BdG Hamiltonian for the dynamics of 
fermionic quasiparticles deep inside the superconducting 



state of a superconductor 
i/ = i(ct, c)U 



S A 

-A* -S^ 



,(A1) 



where His a. AN x AN matrix for a system with N orbitals 
(lattice sites), and c = (c|, C|) is a 2N component vector, 
[c and can be either column or row vector depending 
on the context.] Following the notations in Ref. we 
use two sets of 2 x 2 Pauli matrices to,s,y,z and SQ^x,y,zi 
which act on particle- hole and spin indices, respectively. 
Because of 



^ = (hermiticity) , 
A = — A^ (Fermi statistics). 



the BdG Hamiltonian 
try (PHS) 



(A2) 

1|) satisfies particle-hole symme- 



(«) 



n 



(PHS). 



(A3) 



The presence or absence of TRS and SU(2) spin rotation 
symmetry are represented by 



and 



(6): n^iSyU^i-iSy), (TRS), 



(c) : [H, Ja] = 0, Ja ( "^t 



(A4) 



a = x,y,z, 



-s 
[SU(2) symmetry] 



(A5) 



respectively. 

The ensemble of BdG Hamiltonians ( [All ) with PHS 
condition (a) defines symmetry class D of Altland and 
Zirnbauer. With additional TRS condition (6), the re- 
sulting ensemble of BdG Hamiltonians is called symme- 
try class DHL For both symmetry classes, spin rotation 
symmetry (c) is not necessary. 

The Hamiltonian in symmetry class D can be thought 
of as, because of PHS (a), a single-particle Hamiltonian 
of Majorana fermions. The Majorana structure of the 
BdG Hamiltonians can be revealed by 



where 77 and x ai'c Majorana fermions satisfying 



26,. 



V, [i = 1, 



,2iV), 



(A6) 



etc. 
(A7) 



Then, in this Majorana basis, the BdG Hamiltonian can 
be written as 



where 



X 



If P + S ^i{Q-R) 
2\i{Q + R) P-S 



(A8) 



(A9) 
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and 



P = E-E^ = -P^, 
R = A + A* = -i?^, 



= S + = +Q^, 
S = A- A* = +S'^. 



Then, the AN x iN matrix X satisfies 
A't ^X, X'^ ^ -X. 



(AlO) 



(All) 



These conditions define symmetry class D. On the other 
hand, symmetry class Dili is defined by, in addition. 



iSyX'^ {~isy) ~ X. 



(A12) 



While it is always possible to cast the BdG Hamil- 
tonians into a form of a single particle Hamiltonian of 
Majorana fermions, by rewriting the BdG Hamiltonian 
in terms of the " real" and " imaginary" parts of the elec- 
tron operator, r] and x- there is no natural way in general 
to rewrite Majorana hopping problems as a BdG Hamil- 
tonian. In order to do so, the single particle Majorana 
Hamiltonian must be an even-dimensional matrix, and 
we need to specify a particular way to make a complex 
fermion operator out of two Majorana fermion opera- 
tor. Still, any single-particle Hamiltonian for Majorana 
fermions, with its defining properties (All), can be classi- 



fied in term s of t he presence (class Dili) or absence (class 
D) of TRS ( A12 ) without referring to complex fermions. 

a. off- diagonal block structure of class Dili Hamilto- 
nians Combining class Dili conditions (a) and (6), one 
can see that a member of class Dili anticommutes with 
a unitary matrix t,j;Sy^ 



txSyl~CSytx ■ 



(A13) 



In this sense, class Dili Hamiltonians can be said to have 
a chiral structure. In order to compute the winding num- 
ber V, defined for class DIII Hamiltonians in three spatial 
dimensions, it is necessary to go to a basis in which the 
chiral transformation t^Sy is diagonal. We can find such 



a basis as follows: we first rotate t^ 
by a unitary transformation 



t-, and s„ 



Wi 



1 



1 



V2 



{to + ity) —j= (so - is^) , 



(A14) 



I.e., Wit^Wl = -i^, and W^SyWl = -s^. We then 
exchange the 3rd and 4th entries by a unitary transfor- 
mation W2, 



W-2WltxSyW\wl 



tnS, 



(A15) 



Further exchanging the 2nd and 3rd entries by a unitary 
transformation W3 , the combined unitary transformation 
W = W3W2W1 diagonalizes tx:Sy, 



txSy Wt^SyW'' = t, 



So- 



(A16) 



Under the transformation W PHS and TRS transforma- 
tions are transformed as 

t^so ^ Wt^soW^ = -it^s,, (PHS) 

t^Sy^WtoSyW^ =tyS,, (TRS) (A17) 

respectively. Observe the transformed PHS and TRS 
pick up the same sign under matrix transposition 
[{txSzV = +{txS:,Y and {tyS^Y = -{tySz), respec- 
tively] as the original ones [{txSo)^ = +{txSa) and 
{t^SyY' = ^{toSy)'^, respectively]. In this basis, the 
Hamiltonian takes the block-off diagonal form. 



n 



D 
ZJt 



D = -SzD' sz 



(A18) 



This can be further simplified by a unitary transforma- 
tion 



n 



fit 




oxy 



D 
Lit 



st 

xy 

Sxy 



Sxy-^'^sl-y 





(A19) 



where 



*Z — ~''^S,^yS,^y, S^y ^ ^S^; " S . 



V2 



(A20) 



Introducing 



D' -.^ SxyDsxy = -slyD^sly = -{D'f, (A21) 
we finally arrive at 



D' 
D't 



D' = -D 
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